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Abstract 

Several government/commercial off-the-shelf modu- 
lar software packages have been combined to perform 


plinary optimization (MDO). Results obtained for 
the aeroelastic analysis and aerodynamic optimiza- 
tion cases are consistent with expectations and serve 
as a benchmark for future results that will incorpo- 
rate sensitivity derivatives from the computational 


1 Introduction 

Multidisciplinary design optimization (MDO) meth- 
ods developed in the past decade have not been 
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plications typically involve tens to hundreds of vari- 
ables and hundreds of constraints. In addition, a sin- 
gle high-fidelity aerodynamic or structural analysis of 
an aircraft configuration requires several CPU hours 
on a supercomputer Thus traditional optimization 


an optimal aircraft design. In such a scenario the 
computational expense of applying MDO to realistic 
design problems quickly becomes prohibitive. 


The goal of this research project is to address the 
issue of computational expense in MDO through the 
use of computational fluid dynamics (CFD) solvers 
and finite element (FE) analysis/optimization soft 


FE solvers have recently become available due to 
the implementation of three techniques. These are 
(1) the development of stand-alone sensitivity equa- 
tion solvers for computational aerodynamics such as 
SENSE (see Godfrey and Cliff [1]) (2) the use of au- 
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Halt ka and Gurdal [5]). Both of these methods yield 
sensitivity derivative data at a lower computational 
expense than would be incurred if using traditional 
Unite difference-based techniques to estimate sensi- 


In progress toward the goal of high-fidelity MDO 
for aircraft configurations, an initial capability has 
been developed to perform multidisciplinary analy- 
sis for static aeroelastic problems along with single- 
disciplinary aerodynamic optimization In this re 


sis incorporates the loosely coupled or interfaced 
approach described by Borland [6], where separate 
CFD and FE solvers are employed. The coupling 
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between the two disciplines is accomplished through 
a stand-alone software package which integrates the 
aerodynamic pressures to form aerodynamic loads, 
and then transfers the loads to the FE model via a 
suite of interpolation schemes. While such an ap- 
proach is certainly not new (cf., Tzong et al [7], and 
Chen et al [8]), the uniqueness of this research stems 
from (1) the use of government/commercial off-the- 
shelf software for the CFD solver, FE solver, grid gen- 
erator, and loads transfer, and (2) the goal of using 
available sensitivity derivatives to reduce the compu- 
tational expense of high-fidelity aeroelastic analysis 
and optimization problems. 

This paper contains the background information 
on the suite of software developed for this project 
along with initial results for aeroelastic analysis and 
aerodynamic optimization cases performed without 
incorporating sensitivity derivatives. These data 
form the benchmark results to which future analy- 
sis and optimization cases, performed with sensitivity 
derivatives, will be compared. The methods incor- 
porating sensitivity derivatives in multidisciplinary 
analysis and multidisciplinary optimization are de- 
scribed in Sections 2.2 and 2.3. 

The remainder of this paper is arranged in the fol- 
lowing manner. In Section 2 the aeroelastic analysis 
problem is introduced. Section 3 contains a descrip- 
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presented in Section 5. A summary of this work along 
with items for future work are contained in Sections 
6 and 7, respectively. 

2 Aeroelastic Analysis with 

Sensitivity Derivatives 

2.1 Aeroelastic Problem Formulation 

In static aeroelast.icity the aircraft designer seeks to 
predict the deflections which an elastic structure will 
undergo when subject to time-invariant aerodynamic 
loads. The aerodynamic component of this aero- 
st.ructural interaction is modeled as a solution of the 
governing fluid flow equations, R, in the form 

R[Q(X a ),X a (D)] = 0, (1) 

where D is a set of design parameters describing the 
aircraft geometry, X a is the discretized representa- 
tion of the flow field surrounding the aircraft model, 
and Q is the solution of the state (flow field) vari- 
ables on X a . The vector of aerodynamic forces on 


the surface of the aircraft , F a , is represented as 

F a = F a [Q(X a ), X a (D)]. (2) 

The structural component of the aero-structural in- 
teraction is modeled as 

F s — F s [F a , X S (D)], (3) 

and 

F s = Ku (4) 

In Equation 3, F s is the vector of applied forces (in- 
cluding structural weight) on the discretized struc- 
tural model, X s . In Equation 4, K is the stiffness 
matrix of the discretized undeformed structure, and 
u is the unknown vector of nodal displacements due 
to the applied forces. 

In general, the discretization of the solid bound- 
ary of the aerodynamic grid does not coincide with 
the discretization of the outer surface of the struc- 
tural grid. Thus, a mechanism is needed to transfer 
the aerodynamic loads from X a to X s . In addition, 
another mechanism is needed to map the structural 
deflections, u, to the design parameters, D. Once 
these are available Equations 1-4 may be solved us- 
ing an iterative scheme. 

Note that while nonlinear aerodynamic loads are 
used in this study, the structural model is evaluated 


during the iterative process of solving Equations 1 
- 4. Consequently, during each iteration the vector 
of structural loads, F s , is applied to the undeformed 
structural model. The static wing deflections for a 
HSCT at cruise conditions should exhibit linear be- 
havior (i.e. , small wing tip deflections with respect 
to the wingspan). This assumption of linearity may 
not be valid for off-cruise conditions which are not 
considered in this research. 

2.2 Sensitivity Derivatives in Analysis 

Much of the computational expense in repeatedly 
solving Equations 1-4 occurs in using a CFD 
solver to satisfy Equation 1 as D varies. It is here 
that the sensitivity derivatives are useful in reducing 
the computational expense of high-fidelity aeroelas- 
tic analyses. This is accomplished by approximat- 
ing Q[X a (D + AD)] rather than solving Equation 1 
for D + AD. This approximation, denoted by Q, is 
performed using a first-order Taylor series expansion 
around D of the form 

Q[X a (D + AD]psQ[X a (D)] + ^gAD, (5) 
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where via the chain rule 

dq L= dq L dx^ 

<9D 9X a 3D U 

The sensitivity derivative matrix 3Q/3D may be 
extracted from the output of CFD solver such as 
CFL3D.ADII [2, 9], where CFL3D.ADII is a version 
of CFL3D [10] to which the automatic differentiation 
tool ADIFOR has been applied. In CFL3D.ADII the 
sensitivity derivative matrix is obtained at the same 
time the governing fluid flow equations are solved. 
Alternatively, the sensitivity derivative matrix may 
be obtained using a post-CFD sensitivity equations 
solver such as SENSE [1], In SENSE, the user pro- 
vides the state variables obtained from a CFD solu- 
tion along with a perturbed aerodynamic grid, and 
SENSE returns the sensitivity derivative matrix. 

The approximation of the state variables, Q, is 
then updated as D changes for each iteration of the 
aeroelastic analysis process. Eventually the Taylor 
series approximation of the state variables will de- 
grade and one must solve Equation 1 to obtain accu- 
rate state variables. Thus, a sequence of Taylor series 
approximations of Q is performed until convergence 
is attained. 

Such sequential approximation methods have 
been used by numerous researchers in the past in- 
cluding McQuade et al [11], and Hutchison et al [12]. 
However, these and similar approaches have used fi- 
nite difference techniques to estimate the sensitivity 
derivatives, rather than analytically obtained sensi- 
tivity derivatives. 

2.3 Sensitivity Derivatives in MDO 

Available sensitivity derivatives also permit the use 
of high-fidelity analysis tools in multidisciplinary op- 
timization, in a manner similar to their use in high- 
fidelity analysis described above. One such use of 
the sensitivity derivatives is in the formulation of the 
Global Sensitivity Equations (GSE) developed by So- 
bieski [13]. In a coupled aero-structural problem, the 
aerodynamic derivatives in the GSE matrices would 
be obtained using either CFL3D.ADII or SENSE, and 
the structural derivatives would be obtained from a 
FE solver such as MSC/NASTRAN [14] or GENE- 
SIS [15], both of which calculate sensitivity deriva- 


then be used for optimization studies such as those 
performed by Barthelemy et al [16] and Dovi et al 
[17]. Note that both of these studies employed fi- 


nite difference methods to estimate the terms in the 
GSEs. 

Another MDO formulation that would benefit 
from the use of sensitivity derivatives is the collab- 
orative optimization (CO) strategy of Kroo et al [18]. 
The CO method employs a bilevel optimization tech- 
nique in which a system-level optimizer ensures com- 
patibility of the competing disciplines (e.g., aerody- 
namics, structures, and performance) by prescribing 
target values that each discipline-level optimizer must 
meet. Thus, the discipline-level optimizers must be 
computationally efficient since they are implemented 
numerous times during a CO execution. 

For both the GSE-based and the CO strategies 
described above, the computational expense of using 
finite difference methods and high-fidelity CFD and 
FE solvers is prohibitive. The availability of sensi- 
tivity derivatives has the potential to make the use 
of high-fidelity solvers in the GSE-based and CO ap- 
proaches computationally affordable. 

3 Analysis Models 

3.1 Baseline HSCT Model 

The HSCT model used in this research is a generic 
wing/fuselage configuration extensively studied by 
numerous researchers including Barthelemy et al [16], 
Dovi et al [17], Hutchison et al [12], and Giunta et 
al [19]. The mission performance, and hence the 
economic viability, of such an aircraft hinges on the 
trade-offs among the often competing engineering dis- 
ciplines involved in aircraft design. Thus, the HSCT 
is particularly well-suited to a multidisciplinary de- 
sign approach. In this research, the inter-disciplinary 
couplings are limited to those of aerodynamics and 
structural mechanics. 

The original vehicle used in this study was sized 
to accommodate 252 passengers for a range of 5,500 
nautical miles, and a cruise Mach number of 2.4. [16]. 
Here, the baseline HSCT model is simplified to a 
wing/fuselage configuration. That is, the horizontal 
tail, vertical tail, and engine nacelles are removed. 
While simplified, this model contains sufficient com- 
plexity to be of interest for aeroelastic analysis and 
aerodynamic optimization. Details on the aerody- 
namic and structural models of this HSCT are pro- 


The baseline HSCT configuration is parameterized 
using a technique originally developed by Hutchison 
et al [12] to describe generic aircraft wing and fuselage 
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shapes. For this research project, modifications to 
Hutchison’s technique have been made to accommo- 
date variable camber and twist along the span. In this 
study the area-ruled fuselage and the wing/fuselage 
fairing are held fixed and 64 parameters are used to 
model the variations in planform and wing shape. 

( 


the wing break ( Yb reak ), and the y-location of the 
wing tip (Y tip ). Completing the planform specifica- 
tion is the root chord with a fixed length of 144.87 

ft. 


parameter ( R p ), and the chordwise location of max- 
imum thickness for the entire wing (C ma x. thk .)• hi 
addition, the t/c ratio of the wing root is fixed at 
2.96 percent. In this parameterization the t/c ratio 
f tli i i li 1 b t th t h d 


the radius-to-chord ratio, r/t, given bv 

r/t = 1.1019[(i?p/6)(</c)] 2 , (7) 

where R p is the non-dimensional leading edge radius 
parameter. 

The camber and twist are specified at each of the 
13 airfoils using a cubic polynomial of the form 

Zk = ci t k(x/c) + c 2 ,k(x/c) 2 -I- c 3jk (x/c) 3 , (8) 

where ( x/c ) is the normalized chordwise distance 
specified at 20 locations along each chord, and k = 
1, . . . , 13. Here the 39 polynomial coefficients control 
the spanwise camber and twist distribution 


ment of the wing during an aeroelastic analysis. 

While this 64 parameter model is a simplified rep- 
resentation of a wing/fuselage airframe, this param- 
eterization allows the exploration of a wide range of 


4 Analysis and Optimization 
Software 

4.1 Software Overview 

The software suite developed in this research in- 


sign Optimization) [20] developed at NASA Langley, 
(2) the Euler/Navier-Stokes solver CFL3D [10] also 
developed at NASA Langley, (3) the finite element 
analysis and optimization package GENESIS [15] de- 
1 db VMAE i i I d (4) h 


stitute and funded by the U.S. Air Force. 

This collection of G/COTS software is loosely 
coupled using UNIX shell scripts and some in-house 
software created to parse input and output files as 
d d Th f g t d i ti 


Aeroelastic analysis is performed by coupling the 
software modules in an arrangement depicted by Fig- 
ure 2 where the file interface routines have been omit- 
ted for clarity. The arrows in Figure 2 denote the 
passing of data among the analysis routines. 

In Figure 2 the box labeled ’’Geometry Manipula- 
tor” contains the in-house developed software which 
translates the parametric representation of the HSCT 
model into hie formats for the G/COTS modules. 
For input into CSCMDO an aerodynamic surface 
grid is generated in the HESS format, and for input 
into GENESIS a finite element model is created in 
the MSC/NASTRAN Bulk Data format, with which 
GENESIS is compatible. Ideally, this geometry trans- 
lation would be performed with G/COTS software 
such as Pro-Engineer, AutoCAD, or CATIA. How- 
ever, a number of constraints made it impractical to 


( ) 

and Sun workstations. However, all CFL3D and 
GENESIS analyses are performed exclusively on SGI 
workstations to permit accurate comparisons of ex- 
pended computational resources. Timings reported 
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all calculations performed thus far have been done 
serially. Future efforts will incorporate parallel com- 


Prior to performing an aerodynamic analysis, a vol- 
ume grid is generated based on the aerodynamic sur- 
face grid using the G/COTS module CSCMDO. This 
software tool is an automated, multi-block, struc- 
tured three-dimensional volume grid generator de- 
veloped by members of the Geometry Laboratory 
(GEOLAB) at the NASA Langley Research Cen- 
ter. In addition to grid generation, CSCMDO also 
has the capability to modify an existing volume grid 
based on perturbations to the underlying surface grid. 
It is in this capacity that CSCMDO is used here, 
where a surface grid and a volume grid for the base- 
line HSCT were provided to the authors by GEO- 
LAB. The user then provides a perturbed aerody- 
namic surface grid to CSCMDO which redistributes 
the grid points throughout the computational vol- 
ume using a transfinite interpolation scheme. The 
volume grid generated by CSCMDO is written to a 
file in PLOT3D format, which is a common format 
used by numerous CFD solvers. This efficient grid re- 
generation capability facilitates the use of CSCMDO 
in aeroelastic analysis and aerodynamic optimization 
where the volume grid must be updated repeatedly 
as the structure deforms or is modified during the 
optimization process. 

The CFD solver used in this study is the G/COTS 
program known as CFL3D. This software was de- 
veloped at the NASA Langley Research Center and 
currently is maintained by members of the Aero- 
dynamic and Acoustic Methods Branch. CFL3D 
is a time-dependent, Reynolds- averaged, thin-layer 
Navier-Stokes flow solver for use with two- or three- 
dimensional structured grids. A third-order accu- 
rate finite-volume spatial discretization is employed 
in CFL3D, along with both mesh sequencing and 
multigrid techniques for convergence acceleration. In 
addition, multiblock connection strategies are possi- 


tions are solved on a volume grid with a C-0 dis- 
cretization having the dimensions 121 x 41 x 61, in 
the streamwise, circumferential, and surface normal 
directions, respectively. The volume grid is divided 


is a total of approximately 300,000 grid points in the 
model (see Figure 3). 


conditions correspond to those obtained from the 
U.S. Standard Atmosphere tables listed by Shev- 
ell [24] for an end-of-cruise altitude of 63,000 ft. 
Here, freestream density, pressure and temperature 
are 0.000198 lb ■ sec 2 /ft*, 130.86 lb/ ft 2 , and 392 °R, 
respectively. The reference area used in calculating 
Cl and Cd is based on the planform area of the 
HSCT configuration. For the baseline HSCT half- 
span model the reference area is 4609.65 ft 2 . In this 
model the coordinate axes are aligned such that the 
x-axis points in the positive streamwise direction (i.e. , 
from nose-to-tail) and the y-axis points in the span- 
wise direction from the fuselage centerline to the wing 
tip. 

For all of the calculations presented here the 
angle-of-attack is held constant at 3.5° with zero 
sideslip. Thus, symmetry in the x-z plane permits the 
use of a half-span model, with appropriate bound- 
ary conditions along the plane of symmetry. The 
flow tangency condition is enforced along the solid 
surfaces of the wing/fuselage model and freestream 
conditions are imposed on the outer boundary of the 
computational grid. 

During the initial solution of the Euler equations 
for the aerodynamic model, the entire flow field is ini- 
tialized to freestream conditions. Both mesh sequenc- 
ing and two-level V-cycle multigrid are employed to 
accelerate convergence. The initial flow field solution 
is obtained after 400 iterations which is sufficient to 
converge the residual error parameter in CFL3D by 
five orders of magnitude. The typical computational 
cost for this initial solution is approximately 1.3 CPU 
hours. 

Subsequent solutions of the Euler equations are 
obtained using the restart capability in CFL3D. Here, 
the flow field values are initialized to converged flow 
field values of a previous CFL3D solution instead of 


The convergence history and analysis results of 
CFL3D are output in a variety of files and formats. 
One output option is to generate files in the PLOT3D 
format which contain the grid point locations on the 
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4.3 Structural Model and Analysis 

The FE solver used in this study is GENESIS which 
performs structural analysis and optimization, where 
both element sizing and shape optimization are possi- 
ble. The model definition format in GENESIS is iden- 
tical to that ofMSC/NASTRAN, which is somewhat 
of a de-facto standard in structural analysis. The 
similarity between GENESIS and MSC/NASTRAN 
is useful since many software packages are written to 
read FE models in NASTRAN format 


Balabanov et al [26] ). A copy of this structural model 
generator was provided to the authors and was modi- 
fied for use in this research. The wing/fuselage model 
of the HSCT is comprised of a fixed number and ar- 
rangement of spar and rib elements. The layout of the 
structural elements is based on the planform and air- 
foil thickness parameters used to represent the HSCT 
configuration. 

The finite element model of the HSCT configura- 
tions has 206 nodes and 1130 elements with a total of 
1254 degrees-of- freedom. Note that due to structural 
symmetry only the starboard portion of the model is 
constructed. The FE model contains triangular mem- 
brane elements for the fuselage and wing skins, along 
with rod elements for the spar cap and rib caps. The 
spar and rib webs are modeled with a combination of 
shear panels and rod elements. The material for all 
structural elements is titanium alloy Ti-6A1-4V. The 
structural layout of the baseline HSCT wing/fuselage 
model is shown in Figure 4. 

Some sources of non-structural wing weight are 
incorporated into the half-span FE model. These 
sources include 44,580 lb for engine weight, 7,800 lb 
for the main landing gear weight, and 19,800 lb of 
reserve fuel stored in the wings. These data were ob- 
tained using the FLOPS (Flight Optimization Sys- 
tem) [27] program to analyze the baseline HSCT 
configuration considered here. Note that the engine 
weight was concentrated at four FE nodes on the 


surfaces. However, no attempt was made to remove 
structural elements to model the landing gear bay. 
The fuel weight was equally distributed to 118 FE 
nodes in the inboard wing section and in the leading 
edge st.rake. 

Prior to an aeroelast.ic analysis the spar, rib, and 
skin elements are sized using GENESIS to meet von 


Mises yield stress criteria and local buckling con- 
straints under load conditions approximating a 2.5 g 
pull-up maneuver at Mach 2.4. In this element sizing 
there are 40 variables where 26 of these are skin panel 
thickness parameters (i.e., 13 panels on each of the 
upper and lower body surfaces), 12 variables define 
spar cap areas, and two variables describe the rib cap 
areas. Typically, this sizing optimization requires 15- 
20 CPU minutes on a Silicon Graphics workstation. 

After the initial sizing in GENESIS the element 
di i h Id fi d i th b t 


ESIS structural analysis requires approximately 50 
CPU seconds. 

4.4 Aerodynamic Loads Transfer 

The software used to calculate the aerodynamic loads 
and to transfer them to the structural model is known 
as the Fluids and Structures Interface Toolkit (FA- 
SIT) developed by Smith et al [21, 22]. Currently, 
FASIT is maintained by the Air Force Research Lab- 
oratory at the Wright-Patterson Air Force Base. 

FASIT is ideally suited for the loosely coupled ap- 
proach to aeroelasticity, as it accepts several widely- 
used file formats employed by CFD and FE solvers. 
In this research, the PLOT3D format is used to in- 
put the grid and flow field data from CFL3D and the 
NASTRAN Bulk Data format is used to input the FE 
surface grid. FASIT provides a suite of interpolation 
schemes for transferring the aerodynamic loads from 
the CFD surface grid to the FE surface grid. The thin 
plate spline method of Duchon [28] was used in this 
study, as recommended in the FASIT User’s Manual 
[22]. Note that the interpolation methods in FASIT 
conserve the total force and moments for all three 
axes when transferring loads from the CFD grid to 
the FE grid. 

The loads transferred from FASIT to the FE sur- 
face grid are written in the NASTRAN Bulk Data 
format. Since this NASTRAN format is compatible 


4.5 Structural Deformation 

From GENESIS an output file is obtained which con- 
tains the nodal displacements of the FE surface grid. 
A simple Fortran program was written to transfer the 
nodal displacements to the original (undeformed) FE 


6 

American Institute of Aeronautics and Astronautics 



surface grid to obtain the grid coordinates of the de- 
formed structure. A Mathematica [29] program was 
written to calculate the change in camber/twist vari- 
ables (see Equation 8) and to extract the z-ordinates 
of the leading edge at each of the 13 airfoil sections 
which define the wing. The output from the Math- 
ematica routine is a file containing the parameters 
which are needed by the geometry manipulation soft- 
ware to begin the next iteration in the aeroelastic 
analysis. 

4.6 Aeroelastic Analysis 

Due to the nonlinearity of the aero/structural in- 
teraction the aeroelastic analysis involves an iter- 
ative scheme whereby the aerodynamic and struc- 
tural analyses are performed repeatedly until both 
the aerodynamic loads and the structural deflections 
reach convergence (see Figure 2). Following the 
methods of Chipman et al [30] and Tzong et al [7], 
a constant factor under-relaxation method is used to 
accelerate convergence of this iterative process. A 
value of 0.7 for the relaxation parameter is used for all 
aeroelastic calculations described here. The compu- 
tational cost of the aerodynamic analyses dominates 
the total CPU time involved in obtaining a converged 
aeroelastic analysis. 

4.7 Aerodynamic Optimization 

The modularity of the CSCMDO and CFL3D codes 
permits a simple reorganization of the software for 
use in aerodynamic optimization. Here, the commer- 
cial optimization software package DOT (Design Op- 
timization Tools) [31] is coupled to the aerodynamic 
analysis software using an interface program as shown 
in Figure 5. The sequential quadratic programming 
(SQP) method in DOT is used for the aerodynamic 
optimization case presented here. 

5 Results 

5.1 Aeroelastic Analysis 

The static aeroelastic analysis case is performed to 
calculate the external shape of the deformed HSCT 
structure under 1.0 g Mach 2.4 cruise conditions. Re- 
call that a constant factor under-relaxation method, 
with a value of 0.7, is used to accelerate conver- 
gence of this iterative process. Typically, convergence 
is obtained in six iterations (2.7 CPU hours) when 
under- relaxation is used as compared to 19 iterations 
(6.1 CPU hours) without the relaxation method as is 


shown in Figure 6. The case with relaxation was ter- 
minated when the wing tip deflection was less than 
0.05 ft, while the case without relaxation was termi- 
nated after 19 iterations due to lack of convergence. 

Results for the aeroelastic analysis with relaxation 
are shown in Figures 7 and 8. In Figure 7 the final 
deformed wing is shown in comparison to the solid 
outline of the undeformed wing. Note that the z- 
axis in this figure has been scaled to exaggerate the 
wing deflection for easier viewing. Figure 8 shows 
the aeroelastic deformation of the HSCT wing at the 
wing tip, wing break ( y = 33.7 ft), and wing root. 
At the 1.0 g cruise conditions considered here the z- 
direction wing tip deflection of approximately 1.5 ft 
is reasonable. For this HSCT wing/fuselage model 
at a fixed angle-of-attack of 3.5° (relative to the fuse- 
lage centerline) , this static aeroelastic analysis reveals 
a reduction in lift of about 5.7 percent compared to 
the lift produced by the rigid wing/fuselage model. 

5.2 Aerodynamic Optimization 

The aerodynamic optimization problem utilizes a ten- 
variable subset of the 64 parameters which describe 
the wing planform and shape. In particular, the six 
planform variables and the four thickness variables 
described in Section 3.2 are used in this sample prob- 
lem. Note that the camber parameters are held fixed 
at their nominal values throughout the optimization 
process. The initial, minimum, and maximum values 
for the 10 variables are listed in Table 1. 

The optimization problem is formulated to mini- 
mize the drag coefficient subject to two constraints: 
(1) lift > 210,000 lb, and, (2) a geometric constraint 
on the chord length at the wing leading edge break 
location. The geometric constraint prevents the op- 
timizer from creating a planform shape where the in- 
board trailing edge angle is positive and the outboard 
trailing edge angle is zero or negative. Such a trail- 
ing edge is undesirable from a structural standpoint 
because of the complexity of the load path (i.e. , spar 
layout) that it creates. 

The results of the optimization are listed in Ta- 
ble 2 and demonstrate a reduction of three drag 
counts (A C'd = 0.0003) between the initial and op- 
timal HSCT configurations. This corresponds to a 
decrease of approximately 5.8 percent of the total in- 
viscid drag. 

The computational expense of this optimization 
is about 16 CPU hours on a single-processor SGI 
R10000 workstation where virtually all of the expense 
is in the CFD analyses. The convergence history for 
both C'd and the C'l constraint is shown in Figure 9. 

A planform view of the initial and optimal wing 
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planforms is shown in Figure 10, and the wing tip, 
wing break, and wing root airfoils are shown in Fig- 
ure 1 1 . Recall that camber variables are not included 
in this optimization problem. Therefore, all of the 
wing shape differences are due to differences in the 
10 planform and thickness variables. 

6 Conclusions 

The goal of this work is the development of a high- 
fidelity MDO capability for aircraft configurations. 
Progress toward this goal is shown here through the 
use of modular G/COTS software for static aeroe- 
lastic analysis and aerodynamic optimization cases 


derivatives, obtained from the CFD and FE solvers 
in the software suite, into multidisciplinary analysis 
and optimization cases. 

7 Future Work 

A number of items remain for current and future 
work on this project. The first of these is the use 
of sensitivity derivatives in a sequential approxima- 
tion strategy (described in Section 2.2) for the aeroe- 
lastic analysis and optimization cases. Subsequently, 
the sensitivity derivatives will be incorporated into 
the GSE-based and CO frameworks for MDO as de- 
scribed in Section 2.3. In conjunction with these ac- 
tivities, parallel computation will be used to reduce 
the wall-clock computational expense of these MDO 
cases through the use of the parallel computing ca- 
pabilities in CFL3D and GENESIS. 
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Table 1: Initial values for the 10 HSCT wing variables 
and the allowable range on each variable. 


Variable 

Initial 

Minimum 

Maximum 

A/ 

74.0° 

73.0° 

75.0° 

A o 

45.0° 

40.0° 

50 0° 

Cbreak 

42.4 ft 

40.0 ft 

45.0 ft 

Ctip 

9.3 ft 

7.0 ft 

10.0 ft 

Ybreak 

28.6 ft 

25.0 ft 

30.0 ft 

Ytip 

67.3 ft 

65.0 ft 

70.0 ft 

(t / cjbreak 

2.15% 

1.50% 

3.00% 

/ c)tip 

2.36% 

1.50% 

3.00% 

Rp 

4.5 

3.5 

5.5 

Cmax. thk . 

45.0% 

30.0% 

48.0% 

Table 2: Initial and optimal values for the 10 HSCT 

wing variables. 



Variable 

Initial 

Optimal 

A 

I 

74.0° 

75.0° 

A 

Q 

45.0° 

45.0° 

Cbreak 

42.4 ft 

44.7 ft 

a 

tip 

9.3 ft 

10.0 ft 

Ybreak 

28.6 ft 

28.1 ft 

Ytip 

67.3 ft 

70.0 ft 

(t/c) 

break 

2.15% 

1.50% 

(t / Cjtip 

2.36% 

1.50% 

R p 

4.5 

3.5 

Cmax 

\ thk. 

45.0% 

46.0% 

C D x 10 4 

52.74 

49.70 

Lift 

202,000 lb 1 

210,000 lb 

Planform Area 

4,610 ft 2 

4,745 ft 2 



Figure 1: Planform variables for the HSCT wing. 



Figure 2: The arrangement of software used to per- 
form static aeroelastic analysis. 



Figure 3: A view of the aerodynamic model of the 
HSCT showing the surface grid on the starboard 
wing, the x — z plane of symmetry, and the exit plane. 
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Figure 6: Convergence history of the aeroelastic anal- 
ysis with and without relaxation. 


Figure 4: A view of the structural model of the 

HSCT showing the wing skin elements (port) and the 
rib/spar elements (starboard). 



Figure 5: The arrangement of software used to per- 
form aerodynamic optimization. 



Figure 7: Orthographic view of the deformed wing 
(mesh) and the undeformed wing (solid outline). 
Note the X:Y:Z scaling of 1:1:2 used to show the wing 
deformation. 
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Figure 8: Airfoil sections at the (top-bottom) wing 
tip, wing break, and wing root for the initial un- 
deformed (dashed) and final deformed (solid) HSCT 
wing shapes. 



Y ( feet) 


Figure 10: A planform view of the initial and optimal 
fJSCT configurations. 
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- 0.001 ° 
- 0.002 


Iteration Number 



Figure 9: Optimization history for the 10 variable 
HSCT wing design. 


Figure 11: Airfoil sections at the (top-bottom) wing 
tip, wing break, and wing root for the initial (dashed) 
and optimal (solid) HSCT wing shapes. 
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